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ABSTRACT 

We demonstrate that the radiation emitted by ultrarelativistic electrons in highly nonuniform, small- 
scale magnetic fields is different from synchrotron radiation if the electron's transverse deflections in these 
fields are much smaller than the beaming angle. A quantitative analytical theory of this radiation, which 
we refer to as jitter radiation, is developed. It is shown that the emergent spectrum is determined by 
statistical properties of the magnetic field. The jitter radiation theory is then applied to internal shocks 
of Gamma- Ray Bursts (GRBs). The model of a magnetic field in GRBs (Medvedev & Loeb 1999, ApJ, 
526, 697) is used. The spectral power distribution of radiation produced by the power-law distributed 
electrons with a low-energy cutoff is well described by a sharply broken power-law: P{uj) oc for 
w <^ ujjm and P{u>) oc u>~^p~^'>/^ for w ^ ujjm, where p is the electron power-law index and ujjm is the 
jitter break frequency which is independent of the field strength but depends on the electron density in 
the ejecta, cojm oc y/n, as well as on the shock energetics and kinematics. The total emitted power of 
jitter radiation is, however, equal to that of synchrotron radiation. 

Since large-scale fields may also be present in the ejecta, we construct a two-component, jit- 
ter -f synchrotron spectral model of the prompt 7-ray emission. Quite surprisingly, this model seems 
to be readily capable of explaining several properties of time-resolved spectra of some GRBs, such as (i) 
the violation of the constraint on the low-energy spectral index called the synchrotron "line of death" , 
(ii) the sharp spectral break at the peak frequency, inconsistent with the broad synchrotron bump, (iii) 
the evidence for two spectral sub-components, and (iv) possible existence of emission features called 
"GRB lines" . We believe these facts strongly support both the existence of small-scale magnetic fields 
and the proposed radiation mechanism from GRB shocks. As an example, we use the composite model 
to analyze GRB 910503 which has two spectral peaks. At last, we emphasize that accurate GRB spectra 
may allow precise determination of fireball properties as early as several minutes after the explosion. 

Subject headings: radiation mechanisms: non-thermal — gamma rays: bursts — magnetic fields 



1. INTRODUCTION 

The conventional paradigm of the generation of radia- 
tion by relativistic electrons in magnetic fields is totally 
based on the theory of synchrotron radiation. We demon- 
strate that this theory is invalid if a magnetic field is tan- 
gled on very short spatial scales and we develop a quanti- 
tative theory of radiation in this case. Apparently, the re- 
quired short-scale fields may naturally be present in astro- 
physical shocks. Here we focus on radiation from gamma- 
ray burst (GRB) shocks, for which a detailed theory of the 
formation of magnetic fields has recently been elaborated. 

The relativistic bla st-wave mo del of cosmological GRBs 
(see, e.g., a review by Piran 1999) explains fairly well many 
observational features of this phenomenon, such as the 
rapid variability of 7-ray flux, the prompt optical flash, the 
[ light curvoa and opoctra of delayed aftorglowo, etc.. Thio 



model interprets the prompt ^ - r&y flash as synchrotron 
radiation emitted by Fermi-accelerated electrons in inter- 
nal s hocks propagating in the ejecta (Rees & Meszaros 
1994) and then Lorentz-boosted to the 7-ray band. 'I'he 



afterglows are explained i n a similar way, as the e mission 
from an external shock (Meszaros & Rees 1993) propa- 
gating into the interstellar medium. To achieve the ob- 
served very high luminosities, the magnetic field in the 



GRB shocks must be of nearly equipartition strength, 
eb — B"^ /Snet ~ 1, where ej is the thermal energy density 
of the shocked material. 

Until very recently, t he equipartition assumptio n was 
completely unjustified. Medvedev & Loeb (1999) have 
shown that the relativistic two-stream instability is capa- 
ble of producing magnetic fields with cb ~ 10~^ — 10^^ 
in both internal and external shocks. Observations of af- 
terglow spectra and light curves yield values of eg from 
- 10-1 - 10-2 for GRB 970508 (jwijers fc Galama 1998 



Frail et al. 1999| ; pranot et al. 1999|) to ^ 10" ^ for 
GRB 990123 and GRB 971214 (palama et al. 1999|). Re- 



cent detection of polarization of the optical aftergl ow of 
GRB 990510 ( povino et al. 1999| ; [Wijers et al. 1999D indi- 
cates that the geometry of the magnet ic field in the shock 
is con sistent with the predicti ons of [ Medvedev fc Loel 
T999|) for coUi mated outflows ( [Ghisellini fc Lazzati 1999 ; 



ruzinov 1999|; [Sari 1999|) 



The magnetic field produced in GRB shocks randomly 
fluctuates on a very small scale of roughly the relativistic 
skin depth, which is ~ 10^ cm in internal shocks, for 
instance. On the other hand, the emitting ultrarelativistic 
electrons have much larger Larmor radii. Therefore, the 
electron trajectories are not helical, as they would be in a 
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homogeneous field. Thus, the theory of synchrotron radi- 
ation derived for homogeneous fields is not applicable and 
the spectrum of the emergent radiation is different. Such 
a situation has never been considered in the astrophysical 
literature. 

In this paper we investigate the effect of small-scale mag- 
netic fields on the properties of radiation. We primarily 
focus on internal shocks, for concreteness. We show that 
there are two regimes, depending on the ratio of parti- 
cle's deflection angle and the relativistic beaming angle. 
Which regime is realized depends on the magnetic field 
properties, B and Ab, but is independent of the particle's 
energy. When deflections are large compared to beaming, 
synchrotron radiation is emitted. Otherwise, when parti- 
cle's deflections are small, a new type of radiation — jitter 
radiation — is produced. A quantitative analytical the- 
ory of this radiation is developed in this paper. For the 
power-law distributed electrons with a cutoff, N{j) (x 7"^ 
for 7 > 7niin, where 7 is the particle's Lorentz factor 
and p is the index, the emergent spectrum has the fol- 
lowing properties. First, the spectral power peaks at the 
so-called jitter frequency, tOjm = ^jilmin), which, unlike 
the synchrotron frequency, is independent of the magnetic 
field strength but, instead, depends on the particle den- 
sity in the shock, lu 



Ti. Second, at low frequencies, 
^ 5: ^jmi the spectral power scales as P(w) oc uj^, in con- 
trast to the synchrotron spectrum, for which P{lu) (x lu^/"^. 
The high-frequency portion is, however, determined by the 
electrons, -P(i^) c>c w"'''"^-'/^. Third, the total (i.e., fre- 
quency integrated) powers emitted in the jitter and syn- 
chrotron regimes are identical. 

Since large-scale fields (e.g., due to a magnetized pro- 
genitor, for instance) may also be present in the shocked 
material, we construct a composite, two-component jit- 
ter -|-synchrotron (JS) spectral model of the prompt 7-ray 
emission. We then compare the predictions of this model 
with presently available data collected (mostly) by the 
Burst And Transient Source Experiment (BATSE) on the 
Compton Gamma- Ray Observatory [CGRO). It turns out 
that the JS model is able to naturally explain some prop- 
erties of the GRB spectra which are inconsistent with a 
simple synchrotron shock model. First, almost a half of 
all BATSE bursts violate the so-called synchrotron "line 
of death" prediction, i.e., their low- frequency spectral in- 
dices aje_greate^_tlmi the maximal admissible value of 
1 /3 ( [Preece et al. 1998 ) . The spectra of these bursts are, 
however, well consistent with the steeper uj^-Xecw of jitter 
radiation. Second, the sharp change of the spectral index 
at a peak frequency seen in some bursts is also consis- 
tent with our model. Third, the JS model theoretically 
supports the fact that some GRBs h ave two spectral sub- 
components ( Pendleton et al. 1994[ ). Fourth, there is an 
indication that spectra of so me GRBs ex hibit spectral fea- 
tures — "GRB lines" (see Briggs 1999 for a review and 
some BATSE candidates). It should be noted, however, 
that the results from Ginga and BATSE are somewhat 
controversial. We demonstrate that a line-like spectral fea- 
ture may be associated with a weak jitter component in 
a synchrotron-dominated spectrum. If a future analysis 
shall confirm that "GRB lines" are real features and not 
instrumental (or other) artifacts, then they provide, to- 
gether with the synchrotron component, precise informa- 



tion about properties of cosmological fireballs just a few 
hundred seconds after the explosion. We illustrate this 
on the example of GRB 910503 which has been observed 
by all four instruments on the CGRO and which exhibits 
a second spectral peak at roughly twice the synchrotron 
peak frequency. A simple fit of a spectral shape readily 
yields the value of the magnetic field, ~ 4 x 10^"* in the 
shock, which is in ag reement with results of a completely 
different analysis by Chiang fc Dormer (1999 ). It should 
be noted that a reliable identification/detection of the jit- 
ter spectral features will provide a direct evidence that the 
magnetic fields in GRBs are generated by the two-stream 
instability, since this is the only presently known mech- 
anism capable of producing small-scale, large-amplitude 
fields in shocks. 

At last, we emphasize that the advantage of our JS 
model is that it was not specially designed in order to ex- 
plain peculiarities of the GRB spectra, but solely to study 
the physical effect of small-scale magnetic fields. The phe- 
nomenon considered in this paper is quite general and, 
clearly, relevant not only to the emission from internal 
GRB shocks. A similar mechanism of emission is expected 
to operate in external GRB shocks and, possibly, in more 
conventional supernova shocks and blazar jets. 

The rest of this paper is organized as follows. A qualita- 
tive consideration of the radiation mechanisms in a nonuni- 
form magnetic field is presented in In §|| the struc- 
ture and properties of magnetic fields in GRB internal 
shocks are discussed. In ^ we present a quantitative an- 
alytical theory of jitter radiation. A two-component, jit- 
ter -(-synchrotron spectral model of the prompt 7-ray emis- 
sion is presented in ^ We compare the predictions of the 
model with recent observational results in §|[ Finally, ^ 
is the conclusion. 

2. RADIATION FROM SMALL-SCALE FIELDS: GENERAL 
CONSIDERATION 

Let's consider a plasma at rest threaded by a magnetic 
field. Let's now consider an ultrarelativistic electron with 
the Lorentz factor 7 3> 1 moving in a magnetic field. Be- 
cause of beaming the emerging radiation is concentrated in 
a narrow cone with the opening angle AO ~ I/7 1 in the 
direction of the particle's velocity. In a uniform magnetic 
field the electron moves along a helical trajectory, so that 
the radiation seen by an observer consists of short pulses 
repeated every cyclotron period. The synchrotron spec- 
trum, therefore, consists of a large number of cyclotron 
harmonics, the envelope of w hich is determined by the in - 
verse duration of the pulses (Rybicki & Lightman 1979). 



The spectrum is peaked near the critical synchrotron fre- 
quency ujc = ^j^eB±/meC, where B± = Bcosx and x is 
the particle's pitch angle. 

If the magnetic field is randomly tangled and the corre- 
lation length is less then a Larmor radius of an emitting 
electron, then the electron experiences random deflections 
as it moves through the fleld. Its trajectory is, in general, 
stochastic. This is similar to a coUisional motion of an 
electron in a medium. Bremsstrahlung quanta are emit- 
ted in every collision. Unlike the bremsstrahlung case, 
here "collisions" are due to small-scale inhomogeneities of 
the magnetic field rather than due to electrostatic fields of 
other charged particles. Since the Lorentz force depends 
on particle's velocity, the emergent spectrum will be some- 
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what different from pure bremsstrahlung. There is also an 
ahernative physical interpretation of the process. For an 
ultrar elativistic electron, the me thod of virtual quanta ap- 
plies (Rybicki & Lightman 1979). In the rest frame of the 
electron, the magnetic field inhomogeneity with wavenum- 
ber k is transformed into a transverse pulse of electromag- 
netic radiation with frequency kc. This radiation is then 
Compton scattered by the electron to produce observed 
radiation with frequency ~ j'^kc in the lab frame. 

Keeping this general physical picture in mind, we now 
analyze the problem in more details. Let's consider a 
nonuniform random magnetic field with a typical cor- 
relation scale Xbj the Larmor radius of the electron, 
Pe = jnieC^ /eB± is less or comparable comparable to 
Xb- The emerging spectrum depends on the relation be- 
tween the p article's deflection angle, a, and the beaming 
angle, A6, ( Landau fc Lifshits 1975| ). For ultrarelativis- 
tic particles and small deflection angles, the latter is esti- 
mated as follows. The particle's momentum is p ~ jmeC. 
The change in the perpendicular momentum due to the 
Lorentz force acting on the particle during the transit time 
t ~ Xb/c is p± ~ FlI ~ eB±XB/c. The angle a is then 
a ~ P±/p ^ eB±XB /jrrieC^ . Thus, the ratio of the two 



angles is 



A0 



eB±Xi 

meC^ 



7- 



Xb 

Pe 



(1) 



It is interesting to note that this ratio is independent of 
particle's energy (i.e., of 7) and is determined by the prop- 
erties of the magnetic field only, i.e., by B and Xb- It is 
more convenient, however, to use the wave-vector, fcs, as 
a measure of the magnetic field scale, instead of ^ k^^. 
We now define the deflection-to-beaming ratio as follows. 



7 



There are two limiting cases. 



a 



(2) 




Fig. 1. — Emission from various points along particle's tra- 
jectory; (a) — ct ^ A8, emission from selected parts (bold por- 
tions) of the trajectory is seen by an observer, (b) — a ^ A9, 
emission from the entire trajectory is observed. 

First, 6 ^ a/ Ad ^ 1, i.e., the deflection angle is much 
larger than the beaming angle, see Figure nk. Then, 
an observer sees radiation coming from short segments 
("patches") of the electron's trajectory which are nearly 
parallel to the line of sight (very much like for pure syn- 
chrotron radiation). The magnetic field in every patch is 
almost uniform but it varies from patch to patch. The ra- 
diation is pulsed with a typical duration Tp l/uJc- The 
characteristic frequency of the observed radiation is thus 
Lo ^ LOc- Note that in the time- or ensemble-averaged spec- 
trum, the instantaneous field, B±, which enters ojc, should 



be appropriately averaged, B ~ {B'±}. In this case the 
emergent radiation is completely identical to synchrotron 
radiation from large-scale weakly inhomogeneous magnetic 
fields. 

Second, 6 ^ a /Ad ^ 1, i.e., the deflection angle is 
smaller than the beaming angle, so that the entire elec- 
tron's trajectory is seen by an observer, as shown in Fig- 
ure |llb. The particle moves along the line of sight almost 
straight and experiences high-frequency jittering in the 
perpendicular direction due to the random Lorentz force. 
We therefore refer the emerging radiation to as "jitter" 
radiation. Its spectrum is determined by random acceler- 
ations of the particle. Let's imagine an electron moving 
ultrarelativistically along the line of sight with a constant 
velocity, the transverse accelerations of the electron are 
small. In the laboratory frame, the electron passes through 
the magnetic field inhomogeneities having a typical scale 
Xb ^ k^^ with the velocity c. In the particle's frame 
(i.e., where its parallel velocity vanishes), the field corre- 
lation scale is A'^ ^ Xb/j ^ {kBj)~^ due to the Lorentz 
transformation. The electron's perpendicular acceleration 
changes significantly during r ~ X'^/c {jkBc)~^, so 
that the characteristic frequency of the emitted radiation 
is ujjo ~ ~ ^ksc. In the laboratory frame, this fre- 
quency is boosted to ld-j = ^ojjo- Thus, the spectrum of the 
emergent radiation is peaked at the frequency ujj ~ ^^kBC. 
This frequency is higher than the critical synchrotron fre- 
quency in the uniform magnetic field of the same strength, 
LOc*, namely 

— ~ - <5 « 1, (3) 

as follows from equation (||). 

We should warn here that despite apparent similarity of 
the jitter and free electron laser emission mechanisms, they 
are quite different. The wiggler field in free electron lasers 
is appropriately adjusted for the electron motion to be in 
phase with the produced radiation field to emit coherent 
radiation. Jitter radiation is, in general, incoherent. 

3. THE STRUCTURE OF THE MAGNETIC FIELD IN GRB 
SHOCKS 

To proceed further, a model for a magnetic field in 
GRB shocks is required. We use the only presently avail- 
able quantitative theor y of the magnetic field ge neration 
in shocks proposed by Medvedev fc Loeb (1999| ). To be 



specific, we focus on internal shocks which produce 7-ray 
emission. External shocks which are responsible for the 
delayed afterglows may be treated similarly and will be 
considered in a future publication. 

Shock fronts are shown to be natural sites of the mag- 
netic field generation. Right before a shock, the inflowing 
(in the shock frame) bulk plasma particles meet the out- 
flowing particles which were reflected (scattered) from the 
shock. Such a two-stream motion is kinetically unstable. 
The emergent magnetic field is random with zero mean 
and lies in the plane of the shock front, — perpendicu- 
lar to the shock velocity. In principle, all plasma species 
participate in the instability. We assume the protons and 
electrons to be the only species and discuss their contri- 
butions separately. 

It is important to emphasize that the generated mag- 
netic field fills the entire volume of a shock shell and is not 
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located within a thin layer of order several skin depths 
near the front. There is a gas flow though a shock. Be- 
cause of flux freezing the generated magnetic field is trans- 
ported with the shocked material downstream. This mate- 
rial is replenished with a fresh one where a new magnetic 
filed is thus continuously produced. Since the magnetic 
field is long-lived and does not decay in a dynamical time, 
as indicated by numeric al simulations (see references in 
Medvedev fc Loeb 1999 ), this field will be present in the 
entire ejecta. 

3.1. Fields produced by the electrons and protons 

In this subsection we briefly remi nd main results of the 
theory of Medvedev fc Loeb (1999| ) for future reference. 
Since electrons are light, the instability induced by them is 
rapid: the typical e- folding length (i.e., the e-folding time 
times the shock speed) is much smaller than the charac- 
teristic shock thickness determined by the Larmor radius 
of heavier protons. Therefore, the magnetic field energy 
grows rapidly and reaches the approximate equipartition 
with the electron kinetic energy. 



VelintmeC n ■ 



etrje 



(4) 



i.e., ese = Bl/^iret = {me/mp)rie, where rje is the ef- 
ficiency factor for the electrons which incorporates un- 
certainties due to the nonlinear phase of the instabil- 
ity, one infers from numerical simulations that generically 
rje ~ 0.1 — 0.01, 7int is the relative Lorentz factor of two col- 
liding shells which produce an internal shock, 7int — 
and n is the number density of particles in the expanding 
shell, before the shock. The spatial correlation scale, kse, 
of the field behind the shock is 



„i /4_l/2 



(5) 



where = ATre^n/me is the electron plasma frequency 
squared, 7e is the initial effective thermal Lorentz factor 
of the streaming electrons, and an extra factor of 47i„t is 
due to the relativistic shock compression. 

The protons may generate magnetic fields too. Since 
they are much heavier than electrons, the spatial coher- 
ence length and the e-folding length are comparable to 
the thickness of the coUisionless shock. Therefore, the 
field does not have enough time to grow during the flow 
transit through the shock. It may however grow behind 
the shock if the two-stream motion of the protons persists 
there. If this is the case, then there are two possibilities. 
First, if there is no energy transfer from the protons to the 
electrons or if it is slow compared to the rate of the field 
growth, then the magnetic field energy may be as large as 



2 

-5— = VplintmpC n ■ 
an 



etVp ~ O.let, 



(6) 



so that the resultant field wiU be Bp /Stt ~ B'^/Stt. Al- 
ternatively, if no magnetic field is generated downstream 
the shock, the strength of the field produced by the pro- 
tons may be orders of magnitude lower. Which case real- 
izes may be learned from numerical particle simulations or 
from observations. We keep esp — Bp/8net as a param- 
eter. The characteristic correlation scale of the generated 
magnetic field is 



k 



Bp 



21/4^1/2^ 



(7) 



where cjpp — Ane'^n/mp, 7p ~ 2 is the initial effective ther- 
mal Lorentz factor of the streaming protons. The com- 
pression factor, 47int, is absent because the field is likely 
produced after the compression occurs. 



3.2. The model 

From equations (^), (P, and (^), 
parameter for the electrons. 



we estimate the S- 



1/2 



(8) 



The exact value of (/) = 2~^/'*(7e/7int)^^^ is somewhat un- 
certain: 7e may evolve during the instability from its ini- 
tial value 7e ~ 2 — 3 to 7e ~ 7int '-^ 3 — 4 due to nonlinear 
effects.^ The numerical prefactor may also be affected. 
Therefore we assume that possible values of Se are in the 
range 1 <^ Se 10~^ (given the uncertainty in rj^ from 0.1 
to 0.01) and generically 5e ~ 0.1. The (5-parameter for the 
protons is 



(9) 



provided rjp ^ 0.1. Second, if the energy transfer is fast, 
which may be the case in fields which are in equiparti- 
tion with the electrons or stronger, then the protons may 
efficiently damp their energy into the emitting electrons, 

^Note, the inflowing electrons are cold; they are not Fermi accelerated 
to operate. Otherwise, no magnetic field is produced. 



unless r]p is too small: r]p <^ 10^^ (i.e., esp ^ lO^^eBe)- 

As will be shown below in §^ a spatial spectrum of the 
magnetic field Be with (5e < 1 is required to calculate the 
spectrum of jitter radiation. This distribution of Be over 
scales is difficult to find from the first principles because it 
is determined by fully nonlinear dynamics of the instability 
process. Some constraints may however be drawn. First, 
the two-stream instability produces magnetic fields in a fi- 
nite range of scales, < k < /ccrit.e, where /ccrit.e ~ kse to 
within a numerical factor of order unity. Second, the field 
grows until it becomes strong enough to defiect the parti- 
cles in the transverse direction by _~ 1 radian on a scale 
of the field coherence length, i.e., B~^ (x pe ^b\- If 
now assume that each Fourier harmonic, Bk, is amplified 
independently of others, we obtain: Bk oc k for k < fccrit.e- 
This is the lower limit: the spectrum of the magnetic field 
can only be steeper than linear. 

In reality, all Fourier harmonics are coupled. Thus, 
when at least one harmonic reaches the sub-equipartition 
strength, the streaming electrons are isotropized by ran- 
dom Lorentz forces. This prevents the growth of other 
harmonics. Therefore, the spectrum of the field will be 
steeper than linear and will have a maximum near kse- 
The simplest possible model is a power-law. 



Bk = 



CBk^', ifO<fc<fcBe, 
0, otherwise, 

yet. Note also that 7e cannot be greater than 7iiit for the instability 



(10) 
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where Cb is a normalization constant and > 1 is a 
power-law spectral index of the magnetic field being a 
free parameter here. The constant Cb may be deter- 

mined using Parseval's theorem,"^ \\Ti\h~*aa l_i^i^B\dx — 
i /o°° dki where L is the system size. Taking into ac- 
count that J B'^dx = BiL = B'^cT, where T is the total 
duration of the pulse, we write 



r2 



TT{2^i+l)cTB^k 



-(2/i + l) 
Be 



(11) 



The magnitude of the field, B^, is given by equation 

Strictly speaking, this model may be used for not too 
small fc's: k > fcmin, where fcmin is set by the condition: 
7/(pefcmin) = 1, sce equation (^. The field harmonics 
with k < fcniin are large-scale ones and they contribute to 
synchrotron radiation. Using equation (g), it is easy to 
obtain that kmin = SgkBe- It is also useful to introduce 
the small-scale field component as follows. 



Bl dk, 



Bis 

m 



1 - (j^'^+i 



(12) 



Hereafter, we omit the subscripts by 5; this should not 
cause any confusion. 



B, 




"■Bp l^min "-Be 

Fig. 2. — The model for the magnetic field in GRB shocks. 

To summarize, we assume the following structure of the 
magnetic field sketched in Figure ||. First, there is a mag- 
netic field produced by the electrons, B^, the magnitude of 
which is given by equation (Q). A large fraction of this field 
contributes to the small-scale component, in accordance 
with equation (|l2|). Its distribution over spatial scales is 
described by a power-law with the index /i, see equation 
(p^). The rest, B^ — Bgg, contributes to the large-scale 
component. Second, there is a magnetic field produced by 
the protons, Bp. This field is a large-scale field. Both B^ 
and Bp are random with zero mean. Third, there could be 
an ordered magnetic field left from a magnetized progen- 
itor, B^,. We define the total large-scale magnetic field as 
follows. 



r)2 7^2 



{Bl - B 



Is) 



Bl = Bl 



Bi 



(13) 



4. RADIATION FROM SMALL-SCALE FIELDS: 
QUANTITATIVE THEORY 

In f|2] we qualitatively demonstrated that there are two 
regimes of radiation. An ultrarelativistic electron prop- 
agating through small-scale fields generated only by the 
electrons {5 <C 1) emits jitter radiation. The electron 
moving through larger-scale fields generated only by the 
protons ((5^1) emits synchrotron radiation. The radi- 
ation spectrum from a magnetic field with a broadband 
distribution over scales is neither of the above and must 
be calculated by appropriate scale averaging of a particle 
t raje ctory. However, for a bimodal field distribution of 
§3.2 such that the magnetic energy in B^ harmonics for 
which 5 ~ 1 (i.e., where the transition from a jitter to syn- 
chrotron regime occurs) is small, the separation of scales 
is possible. The resultant radiation will approximately be 
a composition of the jitter and synchrotron spectra.'* For 
clarity and illustration purposes, we derive the radiation 
properties and spectra in both limits separately. A rigor- 
ous treatment of a general case will be given elsewhere. 
All calculations in this section are done in the reference 
frame of the shocked material, unless stated otherwise. 



4.1. Jitter radiation, (5 <C 1 

For sufficiently small (5's such that a <C A^, the velocity 
V of a particle is almost constant whereas its acceleration 
w = V varies with time. Calculating the Fourier compo- 
nent of the electric field using the Lienard-Wiechart (re- 
tarded) potentials, one arrives at the following expression 
for the total energ y emitted per unit solid angle dVL per 



unit frequency duj (Landau fc Lifshits 1975 
Rybicki fc Lightman 1979| , §3.2 



dW = 



27rc3 
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5 77; see also 



(14) 

where "w^i = J we*" * dt is the Fourier component of the 
particle's acceleration, uj' — (1 — n • v/c), and n is the 
unit vector pointing towards the observer. 

First, we can simplify the vector expression in (p^. In- 
deed, in the ultrarelativistic case, the longitudinal com- 
ponent of the acceleration is small compared to the 
transverse component, w\\/w± ^ 1/7^ 1- There- 
fore V and w are approximately perpendicular to each 
other. Second, the dominant contribution to the integral 
over dfi comes from small angles 6 ^ 1/j with respect 
to the particle's velocity. Therefore, we approximately 
write ~ w (1 - v/c + 6*2/2) =i |w (l - v^c^ + 9'^) = 



9^ -I- 7~2). We now can replace integration over the 
solid angle dfl ~ 9d9d(j) with integration over d(j)dLj' /oj 
and integrate equation (|lj) over the azimuthal angle, 0, 
from to 2ti. The spectral energy finally becomes 



dW 
duj 



27rc3 



1 



2u;'274 



duj'. 
(15) 



where Bq is the fraction of the large-scale field which is 
not produced by the electrons. 

^We use the following definition of the Fourier transform: /(tij) = J f{t)e~'"^'^dt. 

^Physically, a particle moves along a helical trajectory about field lines of a large-scale field. This trajectory is slightly perturbed by high fre- 
quency jittering in the (instantaneous) transverse direction due to a small-scale field component. The intensities of the spectral sub-components 
are determined by the magnetic field energy densities at large and small scales. 



6 



The leading term in the above equation is due to high- 
frequency "jittering" of the electron as it moves through 
the random magnetic field. The second and third terms 
in the brackets are corrections due to the angular struc- 
ture of the radiation field convolved with the relativistic 
beaming. 

The acceleration w is found from the equation of mo- 
tion, p = (e/c)v X B. In general, B may vary both in 
amplitude and in direction. In relativistic GRB shocks, 
radiation is beamed; only a conical section bound by the 
opening angle Ob ~ l/7sh is seen by an observer, where 
7sh is the bulk Lorentz factor of the expanding shell. At 
a particular observing time, mo st of the radiation is com 



ing fr om the b rightened limb (3ari 1998; Panaitescu & 
Meszg xos 1998|). Becau se of relativistic aberration (Ry 



bicki Lightman 1979| ), the shock front is seen almost 
edge-on at the limb. On the other hand, the magnetic 
fields in the shock shell are random but always lie in the 
plane of the shock front (Medvedev & Loeb 199£). Choos- 
ing a coordinate system with x-direction pointing towards 
the observer (in the shock frame) and x-y-plane being the 
shock front plane at the limb, we have v = ci, Bj^ = _By, 
and p ~ "fmeWZ. Then, the Fourier component of w is 



eBk' 



(16) 



For the magnetic field distribution, we use the model (0). 

At last, the spectral power P{uj) = dW/dtdu! (i.e., the 
energy emitted per unit time per unit frequency interval, 
measured in ergs Hz"'^) is obtained dividing the en- 
ergy spectrum dW /duj by the pulse duration, T dRybickj 
& Lig jhtman 1979| ). Combining equations (|l|), ([u]) 
together and substituting into (P_5|), we obtain 



Pico) 



r2 
2uji 



(17) 



where uj/ujj < 2, and r^, — /rrieC^ is the classical electron 
radius. As is expected, the characteristic frequency of the 
emergent radiation is 

LU, = -f'^kBeC = 2^/S'7int7e"'^''^pe, (18) 

cf., equations (||) and (||). The function J is defined as 



J(0 = (2/i+l)e'^ 



/ min 



^' S 



no 



(19) 



the more peaked is the radiation spectrum and the closer 
the peak frequency to ^luj. Note, the discontinuity of the 
slope is artificial. This is due to the fcmin-cutoff because the 
model we use does not continuously interpolate between 
the jitter and synchrotron radiation limits (i.e., between 
6-^1 and 5^1 cases). This discontinuity is less promi- 
nent for large ^'s and small S's. 




Fig. 3. — The arbitrarily normalized power spectra, P{lo), of 
jitter radiation (linear-linear plot) emitted by a single electron 
for several values of the magnetic field spectral index, fi = I 
(dotted), fi = 3 (dashed), and fi = 10 (solid). Here 5 = 0.6 
and we keep the value of the small-scale magnetic field fixed, 
B|s/87r = B^JSn. 

The total emitted power may be obtained by inte- 
grating (|l^) over frequencies. A simpler way is to ex- 
press it in terms of electr on's acceleration, i.e., dW/ dt = 



{2e^-/'^/3c^){wl +-f^wj) ( [Rybicki fc Lightman 19"79| ) and 



substitute w^^ = and Wj_ given by (iq). The result is 



dW/dt 



(21) 



Note, this expression is identical to that for a synchrotron 
radiation in which a uniform field, B^, is replaced with its 
average, B^g = {Bl). 

We are now able to calculate the jitter radiation spec- 
trum from an ensemble of electrons. We assume that elec- 
trons are accelerated in the shock^ to a power-law distri- 
bution N{'j) = Cn"/~^ (where Cn is a normalization con- 
stant) with a minimum Lorentz factor, i.e., 7min < 7 < oo. 
The index p must be p > 2 so that energy does not diverge 
at large 7's. We assume the standard value, ^ p = 2.5 
( Sari, Narayan, fc Piran 1996| ), unless stated otherwise. 



'he spectrum is found from 



where / is the integral,^ /(^) = J (l - C + ^C^) d^, 
and min[a;61 denotes the smallest of a and b. From equa- 
tions (H), (y), (18), we have ujj = ^'^B^./Sm^c. Equation 
(|T^ ) may now be cast into the form. 



Bis 

m 



(20) 



-Pens(w) 



iV(7)P(c^)d7, 



(22) 



noting that ujj = ^^^^(7) oc 7^, i.e., in exact analogy with 
synchrotron radiation, ujc oc 7^ 



27^7! 
' min 



We obtain 







(23) 



This spectrum is shown in Figure |^ for several values of 
/z. In general, the steeper the field distribution, Bf^ cx fc'^, 

^We keep this integral in a general form because it contains logarithmic terms for = 0.5, 1, 1.5. 

® Small-scale magnotY fjclHc prnccnt in tVip ghr.fV provide effective collisions in the otherwise coUisionless plasma and make the Fermi 
acceleration to operate (Medvedev & Loeb 199£). 

^Resent studies indicate that the high-frequency spectral power index in the prompt GRBs is ^ —1 which translates to p 3. We keep 
p = 2.5 for illustrative purposes. 
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where we introduced the characteristic frequency of radi 
ation ujjjn = (^j(7min), cf. equation (1^ 



(24) 



At low and high frequencies, the spectra may be obtained 
analytically as follows, 



Pcns(w) OC < 



-(p-l)/2 



(25) 



At high frequencies, the spectrum is analogous to the syn- 
chrotron case. The electron power-law index determines 
the slope. At low frequencies, the spectrum is linear in 
frequency, in contrast to the synchrotron spectrum which 
is softer and scales as Psyn(w) oc uj^^^. The spectrum of 
radiation from single-speed electrons is peaked at ~ 2LUj 
for large /x's, i.e., for a peaked magnetic field distribution, 
as discussed above (see also Figure^). Therefore, the spec- 
tral break due to the 7niin-cutofF occurs at 



^j,brcak 



2UJ^ 



.,„/2^ = (27/V^)7Ln7int% '/'c^pe 

6.0X 109 7Ln7int7e"'/'%o' Hz, 



(26) 



where riio = n/lO^" cm~'^. The above frequency is calcu- 
lated in the frame of the relativistic expanding shell. This 
frequency is boosted in observer's frame by a factor of 7sh. 

The model of jitter radiation contains two extra param- 
eters, compared to the model of synchrotron radiation. 
These parameters are the deflection-to-beaming ratio, S, 
and the magnetic field index, which determines the 
peakedness of the magnetic field over spatial scales. We 
now investigate how the spectrum of the emergent radia- 
tion depends on these parameters. 
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Fig. 4. — Spectral power P(aj) of radiation (log-log plot, arbi- 
trary units) emitted by the power-law distributed electrons vs. 
^I^cm for (5 = 0.9, 0.7, 0.5, 0.3, 0.1 (from top to bottom). Here 
/i = 1 and p = 2.5. The synchrotron spectrum with the same 
magnetic field strength is shown (dashed curve) for comparison. 
The value of Bg^/Bg is kept fixed for all (5's. 
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Fig. 5. — Spectral power P((i;) of radiation (log-log plot, arbi- 
trary units) emitted by the power-law distributed electrons vs. 
^ I'-^jm for /i = 0, 1, 3, 10 (from top to bottom). Here (5 = 0.7 
and p = 2.5. The synchrotron spectrum with the same magnetic 
field strength is shown (dashed curve) for comparison. The value 
is kept fixed for all ^'s. 



of 



Figure ^ re pre sents the spectral power, P(w), given 
by equations (pO|), (^), as a function of w/wcm, where 
iOcm = '^c(7min), for various values of 5. The dashed curve 
represents the synchrotron spectrum which corresponds to 
the same magnetic field strength. Note, the frequency is 
normalized by the synchrotron break frequency, LOcm (not 
the jitter break frequency, Wjm), to emphasize that, for the 
fixed field magnitude, Pgg, the jitter frequency increases 
with decreasing ^, in accordance with equation (^). One 
can see that jitter radiation is well described by a broken 
power-law given by equation (|2^). Notice also little change 
in the spectrum shape near the break frequency and the 
overall decrease of P(a;) as 5 decreases. 

Figure |^ shows P(lS) vs. uj/u>jrn for a few values of 
the magnetic field index, ^. In general, the ratio Sgg/B^ 
depends on ^. To highlight the effect of the magnetic spec- 
trum shape alone, we keep Bgg/B"^ fixed for all /i's. One 
can see from Figure ^ that has little effect on the overall 
shape of the radiation spectrum. However, as one goes 
from a flat (/i = 0) to a peaked (/i — 10) magnetic field 
distribution, the change in slope at the break frequency 
becomes more abrupt. 

4.2. Synchrotron radiation, (5 3> 1 

For completeness, we also consider the case of (5 ^ 1. In 
this case, a/A9 ^ 1 while still Xb/pl ^ 1 for energetic 
electrons. The spectral power emitte d by a single electron 
IS given as ( [Landau fc Lifshits 1975| , §77 and §74) 



dW 
dt dio 



where Qc 



V37rc72 
i7^e-B-L/™eC, (. 



F 



(27) 



.) denotes the average over 
the electron's trajectory, and "tilde" denotes the local (in- 
stantaneous) quantity. Here F{£_) = ^ K^/^{^') d^' and 
K^/^{x) is the modified Bessel function of 5/3 order. As 
for the spectrum from a plasma, the observed radiation 
conies from regions with different field strength. There- 
fore, one has to replace B± in (G^) with the ensemble av- 
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erage, B 



which rigorously yields the standard 



synchrotron spectrum. 



5. JITTER+SYNCHROTRON MODEL OF GRB EMISSION 

We now use the model of magnetic fields discussed in 
S 3.2 to construct a GRB spectral model All calculations 
are performed in the frame of the expanding shell. The 
transition to observer's frame is obvious. The magnetic 
field in the shock shell was shown to be sub-divided into 
small-scale and large-scale components. The small-scale 
field, for which 6 <^ 1, yields jitter radiation. Its spec- 
tral power is given by equations (J^), (^). The large- 
scale field, for which (5 3> 1, yields s ynchrotron radiation. 
The r elevant theory may be found in Rybicki fc Lightmarj 
19791 ). The composite spectrum, thus, contains both jit- 
ter and synchrotron components and may be schematically 
represented as follows. 



Pj+s{uj) = Pj (w; Bss, 5, + Ps (^; B 



LS 



(28) 



This approximation is valid for the field distribution from 
§ ^.2| unless (5 '--^ 1, as discussed in the beginning of ^ It 
is convenient to normalize frequencies onto the jitter fre- 
which is independent of the magnetic field 



quency, ujj„, 
strength, cf., equation 
quency ratio is then 



(|26|). The synchrotron-to-jitter fre- 



3 Bls 
2 Bss 



(29) 



as follows from equations (||), (p^). 





0.4 




0.2 









-0.2 










O 
J 


-0.4 




-0.6 




-0.8 



-1.5 



-0.5 
Log (o/cOj 



0.5 




Log (B/C0j,„ 

Fig. 7. — A typical composite power spectrum, Pj_|_g(a)), (log- 
log plot, arbitrary units) for the power-law distributed electrons 
with p = 2.5 for /i = 1, (5 = 0.3, and B^/Blg = is shown. The 
synchrotron (long-dashed curve) and jitter (short-dashed curve) 
sub-components are also shown. 

Figure ^ represents the spectrum emitted by single- 
speed electrons for the same values of /i as in Figure 
One can clearly see a sharp feature on top of the broa^ 
synchrotron spectrum. Integrating this spectrum over the 
power-law distribution of 7's, we obtain the composite 
"jitter -f synchrotron" (JS) model of GRBs. A typical ex- 
ample is shown in Figure ^. In general, there are two 
bumps: the sharp one is near the jitter frequency and 
the other, broad bump, is associated with synchrotron 
emission. Depending on the ratio Bj^s/^SS' these bumps 
may overlap to produce either featureless broad or sharply 
peaked spectra, as well. The high-frequency tail always 
scales as Pj+s(a;) oc 0;"^^"^^/^. The dependence of the 
spectrum from S is displayed in Figure M. As an exam- 
ple, we take Bq/Btt = Bgg/Sn and fx = 10 (for such /x's 
BIs/8tt ~ BI/Stt, and, therefore BIs/8tt ~ B^/8tt). A 
sharp jitter feature is easily seen in the spectrum. As S 
decreases, the synchrotron bump moves towards lower fre- 
quencies and decreases in amplitude. The jitter peak de- 
creases even faster and the spectral feature becomes less 
prominent. The other parameter, /i, determines the ratio 
Bj^g/Bgg but, besides this, its effect onto the spectrum is 
weak (cf.. Figure]^) and is, therefore, not shown. 



Fig. 6. — The arbitrarily normalized, composite power spec- 
tra, Pj^g{ijj), (log- log plot) emitted by a single electron for sev- 
eral values of the magnetic field spectral index, A* = 1 (dotted), 
/i = 3 (dashed), and fi = 10 (solid) for S = 0.6 and B'^/B^g = 3. 
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Fig. 8. — The composite power spectra, Pj-|.s(lj), (log-log 
plot, arbitrary units) for the power-law distributed electrons 
with p = 2.5 for 5 = 0.9, 0.3, 0.1, 0.033, -Bg/S" 
M = 10. 



|g = 1 and 



The spectrum of radiation depends on the relative mag- 
nitudes of the large- and small-scale field components. 
Figure ^ represents spectra for /x = 10 and various val- 
ues of the ratio Bg/Egg ~ ^Is/^ss- Figure ||a shows 
the spectral power vs. frequency. For strong large-scale 
fields, the synchrotron spectral component dominates. As 
B\g/Sgg decreases, the synchrotron peak moves towards 
lower frequencies. The amplitude of the synchrotron 
peak decreases too. The position and the amplitude of 
the jitter peak remain almost unchanged. Thus, in the 
limit B^g/Bgg —^ 0, the spectrum becomes purely jit- 
ter. It is illustrative to depict the spectral flux, F{lu) = 
P{uj)/uj [observationally, this quantity is proportional to 
the photon spectral flux, N(E), measured in units: pho- 
tons s~^ cm~^ keV~^], as presented in Figure |b. It IS seen 
that the slope of the fluj^ below the jitter break continu- 
ously decreases as Bj^s/^ss decreases. The synchrotron 
asymptotic slope oc w"^/^ is shown for comparison. It is 
interesting that for Bj^s/^ss ^ 10~^ (the actual value 
depends on 6; the larger 6, the larger the ratio), there is 
a portion in the spectrum which has the power-law index 
being less than —2/3 and approaching zero in the limit of 
the vanishing large-scale field. We discuss this property 
in §^ in the context of the "line of death" for synchrotron 
radiation in GRBs. 

Finally, it is worthwhile to compare the peak spec- 
tral fluxes of jitter, Fj^max = Fj(uijm), and synchrotron, 
-Fs.max = Fsiujcm), radiation. 



F 



J, max 



(30) 



where f{p,fi) is a function of two power-law indices, p 
and fi. The above equation may be readily obtained from 
equations (|3|), (^), and, for instance, (|27|). For /i 1, 
the function / depends on fi only weakly while, given the 
spectrum, p is found by fitting the large frequency slope. 
Therefore, for a given spectrum with two sub-components, 
both significant parameters, namely B^s/^ss Siiid S, are 
uniquely found from equations (p9|) and (|3C|). 
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Fig. 9. — The composite spectra (log-log plots, arbitrary 
unitsj for the power-law distributed electrons with p = 2.5 for 
^o/^SS = 10 - 10"". ^ = 0.2, and ^l = 10. (a) — Spectral 
power Pj+s(tj); (b) — Spectral flux Fj+s(a;) = Pj^si^)/^- 



6. COMPARISON WITH OBSERVATIONS 

In §H we constructed a composite, jitter -|-synchrotron 
model of radiation emitted from internal shocks of 
GRBs, assuming that the magnetic fields are produced 
i n these shocks via the re lativistic two-stream instability 
(Medvedev fc Loeb 199£). Several observational predic- 



tions can now be drawn. We compare them with presently 
available observational data. For future reference, we de- 
fine the photon spectral index, s, as follows. 



F(w) cx to'' 



(31) 



Here we remind that F{uj) = P{lo)/u) cx N{E), where 
N{E) is the number of observed photons per unit time, 
per unit energy range, per unit area. 

We have shown that prompt ORB spectra consist of 
two spectral sub-components, namely synchrotron and jit- 
ter. The synchrotron component is, as usual, well approx- 
imated by the smoothly broken power-law , also referred 
to as the Band function ( Band et al. 199^ ) or the GRB 
function. The jitter component is better approximated by 
a sharply broken power-law with the hard and soft photon 
indexes being equal to s = —(j)+l)/2 and s = respec- 
tively. The position of the jitter peak is independent of 
the magnetic field strength, in contrast to the synchrotron 
peak, but depends on the particle (electron) density in the 
relativistic expanding shell, see equation (26). 



6.1. Low-frequency spectra and the "Line of Death" 



10 



The optically thin synchrotron model of GRBs makes et al. (1998 ) are above t he s — line by mo re than Icr 



a solid prediction. Namely, the soft photon spectral in- 
dex must be in the range —3/2 < s < —2/3, depending 
on the strength of the magnetic field if strong synchrotron 
cooling of the emitting e l ectrons in the fireball is take n 
into account (Katz 1994; 3ari, Narayan, fc Piran 19961 ). 
Thus, the photon index in this model can newer be greater 
than —2/3, creating a t estable "lin e of death" for the syn- 
chrotron shock model ( Katz 19^ ) . There is, however, a 
growing observ ational evidence that many bursts violate 
this prediction (Crider et al. 1997|; ^trohmayer et al. 1998 



Preece et al. 1998 



prontera et al. 1999] ).' For instance 

Preece et al. (1998) have studied time-resolved spectra of 
the bursts collected by the LAD from the BATSE instru- 
ment. They found 23 out of 137 b ursts which violate the 



optically thin synchrotron model. Frontera et al. (1999) 
demonstrated that about 50% of time-resolved spectra the 
bursts observed by the Wide Field Camera on board of 
the BeppoSAX before May 1998 also violate this model. 
The proposed explanations, whic h include Compto n up- 
scattering of low-energy photons ( Liang et al. 1997 ) , syn- 
chrotron self- absorption ( Papathanassiou 1999), and the 
influenc e of the pair annihilation in the fireball photo- 
sphere (Eichler & Levinson 1999), though possible, seem 
rather ad hoc and suffer from drawbacks. The Compton 
up-scattering model strictly requires a single up-scattering 
event per photon which requires the column density to self- 
adjust to a few g/cm^. The self-absorption model results in 
large optical depths and, thus, weak emitted flux and very 
low radiation efficiency. The photospheric model requires 
an extremely lo w baryonic load. For more discussion see 
Meszaros (1999D. 
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Fig. 10. — The power-law index, s, vs. log{Lu/ujjm) for syn- 
chrotron and jitter radiation with fi = 10, <5 = 0.2, B'^^/ B'^g = 
1. The synchrotron and jitter "Unes of death" (LoD) are shown. 



The JS model provides a natural resolution of this puz- 
zle. Figure |l^ represents the photon index as a function 
of a logarithm of frequency for the synchrotron and jitter 
cases* When jitter radiation dominates, the low-energy 
photon index tends towards its asymptotic value of s = 0, 
thus shifting the "line of death" to harder spectra. In prin- 
ciple, the indexes as large as s = 1/2 are allowed by our 
model. Positive s, however, are quite unlikely. These pre- 
dictions are in good agreement with observationa l data. 
Namely, only two bursts out of those studied by Preece 



and all bursts studied by Frontera et al. (1999) except the 
peculiar one, GRB 970111, fall below this line. At last, 
these three bursts are still consistent with s = 1/2 within 
statistical uncertainties. 

6.2. Sharply broken power-law spectra of GRBs 

Most of the observed burs ts are well-fitted by the GRB 
function (Band et al. 1993). However, there are bursts 
with a sharp curvature of the spectrum at the break en- 
ergy. For such bursts a broken power-law model generates 



better fits than the GRB function (see, e.g., Preece et 



al. 1998). In terms of the JS model, there are two cases 
in which such spectra occur. They can be observationally 
distinguished by the value of the soft (i.e., low-energy) pho- 
ton index. First, it is the case of purely jitter radiation: 
B\g/B^g and there is no synchrotron peak (or this 
peak is too weak and is at low energies, outside the detec- 
tor range), see Figure p[ These spectra are flat, s ^ Q, at 
low energies. Second, for large 6 1 and the equipartition 
between the large- and small-scale fields, Bj^s/^ss ~ 
similar spectral shapes are also obtained, see Figure ^. In 
this case, the soft photon index is s ^ —2/3, or it is in the 
range —3/2 < s < —2/3 if s ynchrotron cooling of the clcc - 
trons is taken into account ( ^ari, Narayan, fc Piran 1996 ). 
Thus, we expect fewer broken power-law bursts which have 
the soft photon indexes in the range —2/3 < s < 0. Fig- 
ure 2 of Preece et al. (1998 ) indeed reveals a gap between 
the power-law bursts with s ~ and s < —2/3. How- 
ever, a larger number of such bursts is required to draw 
a statistically significant conclusion. We also should point 
out that if no ordered field is present in the ejecta, then 
Bls ^ Bss is equivalent to B'^/Sn - B^/Stt. This, in 
turn, naturally implies strong energy coupling and rough 
equipartition bet ween the protons and the electrons, as 
discussed in St5.ll 



6.3. Two-component spectral structure of GRB emission 



The analysis done by Pendleton et al. (1994) of the 



bursts from the first BATSE catalog collected by the Large 
Area Detector (LAD) demonstrates that there is a large 
number of bursts which have the high-energy photon in- 
dexes (50 — 300 keV range) being much larger than the 
low-energy ones (20 — 100 keV range). No such behavior 
is observed at other energy ranges. Thus, this result may 
indicate the presence of more than one spectral compo- 
nent in the energy range 40 — 100 keV. (The absence of 
similar behavior at higher energies seems to rule out the 
inverse-Compton origin of a second spectral component.) 
On the other hand, such spectra are completely consis- 
tent with the JS model. They generally correspond to the 
well-separated synchrotron and jitter peaks and are char- 
acterized by small (5's, 5 ^ 0.1, and weak large-scale fields, 
^Ls/^ss ^ 0.01, as is clearly seen from Figure ^ 

6.4. "Lines" in GRB spectra 

Whether emission and/or absorption features, often re- 
ferred to as lines, in GRB spectra are real remains an 
open quest ion for a long time (see, e.g., Briggs 1999 and 
Ryde 1999 for discussion). These features have been ob- 
served, for instance, by the KONUS experiment, Venera 



'Little spikes on the lower curve are numerical artifacts. 
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Fig. 11.— A visual fit of a spectral shape for GRB 910503. The inferred values are: B^/B^g = 7, 5 = 0.07, = 10, p = 3.9. The 
synchrotron (long-dashed curve) and jitter (short-dashed curve) spectral sub-components are also shown. 



mission ( Mazets et al. 1981 ) and by Ginga ( Murakami et 
al. 19|8|). At that time, it was widely believed that these 
features are cyclotron lines due to a highly magnetized 
progenitor. The BA TSE Spectroscopy D etectors (SD) re- 
fined these results. Palmer et al. (1994 ) studied almost 
two hundred bursts detected by SDs. No convinci ng line 
features were found. In a more recent analysis by Briggs 



(19991) of more than a hundred bursts observed by SDs 
before May 31, 1996, about 10 highly significant line fea- 
tures have been found. These are low-energy (~ 50 keV) 
emission features. Clearly, this controversy requires fur- 
ther studies. However, it seems likely that if lines exist, 
they are rare. 

Not too surprisingly, the two-component JS model is 
able to produce spectra with emission line-like features. 
These are not lines in a strict sense because at energies 
higher than the "line" peak, the spectrum has no curva- 
ture and falls down as a power-law with s = — (p + l)/2. 
Few such spectra are seen in Figure ^. These spectra are 
obtained (i) for rather low values oi S {5 <^ 0.3), otherwise 
the synchrotron and jitter peaks overlap, and (ii) for a nar- 
row range of B^s/^SS (^-l ^ ^ls/^SS ^ otherwise 
either the spectral components are too well separated if the 
value of this ratio is small, or the jitter peak is too week 
and unobservable if the ratio is large. We thus conclude 
that the spectra with emission features are not very com- 
mon. They require some parameter tuning and, therefore, 
must be rare. Moreover, inhomogeneities and highly vari- 
able conditions in the GRB shocks may smear out weak 
spectral features completely. 

At last, we present an illustrative example. Figure ^ 
shows the spectrum of GRB 910503 obtained using all ca - 
pability of CGRO^s four experiments ( ^chaefer et al. 19"98|) . 
Here the energy flux vF^, oc E'^N{E) (x u!P{u!) in units 
of (photons cm-2 s'^ keV'^) x (£;/10G keV)^ is plotted 
vs. frequency in units of keV. A sharp second peak at 
^ 2 MeV is clearly seen. The solid curve in the figure 
represents a visual (i.e., not x^) spectral shape fit using 
the JS model. Here we fit the separation of the two peaks 
and their relative amplitudes; we didn't fit the absolute 
flux and position of the synchrotron spectral break which 
depend on the energetics of the fireball. Note that the in- 



tensity of a second, high-frequency (jitter) component is 
low. The values inferred from this fit are the following: 
Big/ Big = 7,5 = 0.07, p = 3.9, and we took = 10. We 
calculate the efficiency of the magnetic field generation in 
the shock from equation (^) assuming = 1/4 as a typical 
value. We obtain r^e ~ 0.08. The total magnetic field en- 
ergy is calculated as eB_= eBe{^ + Blg/Bgg), where ese is 
found from equation (j^ and we used Bss — B^ for large 
/i. We obtain eg — 4 x 10~^, which is in agreeme nt with 
the conclusion drawn by Chiang fc Dermcr (1999[ ) from a 
completely different analysis that the magnetic field in this 
burst is well below the equipartition, es <, 10^^. Their ar- 
gument goes from the fact that the low-energy spectrum 
of this burst is consistent with F{uj) cx w"^/'^. Evidently, 
synchrotron losses are small and the emitting electrons do 
not form a cooled distribution, which otherwise would re- 
sult in F{u}) cx uj~^/'^. This is possible only for very small 
values of eg. 

7. CONCLUSION 

In this paper we have shown that radiation produced 
by relativistic electrons in magnetic fields may differ quite 
substantially from synchrotron. This radiation, referred 
here to as jitter radiation, is produced in the magnetic 
fields which are highly inhomogeneous on very small spa- 
tial scales. Such fields are likely to be present in GRB 
shocks. We developed a quantitative theory of jitter radia- 
tion. Jitter radiation has a different spectrum and its peak 
frequency is independent of the field strength. However, 
the total (i.e., frequency integrated) emitted power of jitter 
radiation depends on the field strength and is exactly iden- 
tical to that of synchrotron radiation. We also constructed 
a composite, two-component jitter -|-synchrotron spectral 
model of the prompt GRB emission. Predictions of this 
model seem to be in excellent agreement with presently 
available data and, likely, resolve some puzzling spectral 
properties of the prompt 7-ray emission. All this, we 
think, strongly supports that (i) the proposed jitter ra- 
diation mechanism operates in astrophysical objects and 
(ii) the magnetic field is generated in shocks by the two- 
stream instability. We emphasize that a reliable identifi- 
cation/detection of the jitter spectral features will provide 
a direct evidence that the magnetic field in GRBs is due to 
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the two-stream instability, since we presently unaware of 
any other mechanism which is capable of producing the 
required small-scale, large-amplitude fields. In general, 
the detection of both spectral components in GRB spec- 
tra would be a powerful and presize tool to investigate the 
properties of cosmological fireballs. 

It is important to emphasize that the phenomenon of jit- 
ter radiation is intrinsic not only to internal shocks. Simi- 
lar conditions (i.e., strong, small-scale fields) are expected 
to occur in external shocks which produce delayed after- 
glows, as well as in more conventional supernova shocks 
and relativistic jets. More speculatively, the above mecha- 
nism could allow to study the magnetic and electric fields 
in reconnection regions (remember, reconnection occurs on 
the electron skin depth scales) and small-scale structure of 
magnetic turbulence and cascade (e.g., in the interstellar 
medium) . 

The theory presented in this paper is only a "first 
step". It is incomplete in a sence that it does not in- 
clude the effects of both ordered and large-scale random 



magnetic fields in a self-consistent way and, hence, does 
not smoothly interpolate between the two extremes of syn- 
chrotron and jitter mechanisms. The synchrotron self- 
absorption has been omitted. The model also does not 
consider some related effects. For instance, it is likely that 
random electric fields (i.e., strong Langmuir turbulence) 
are present in the cosmological coUisionless shocks. These 
fields will definitely affects the radiation spectrum via a 
similar mechanism. It is also unclear now whether par- 
ticle's motion is ergodic in random fields (i.e., whether 
a particle "samples" strong and weak fields statistically 
homogeneously), and what could be the effect of noner- 
godicity on the observed spectrum. All these issues will 
be addressed in future publications. 

The author is grateful to Ramesh Narayan and Norm 
Murray for their interest in this work, various insightful 
comments, and useful discussions and to George Rybicki 
and Avi Loeb for discussions. This work has been sup- 
ported by NSF grant PHY 9507695. 
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